Quantum simulations and experiments on Rabi oscillations 
of spin qubits: intrinsic vs extrinsic damping* 
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Electron Paramagnetic Resonance experiments show that the decay of Rabi oscillations of en- 
sembles of spin qubits depends noticeably on the microwave power and more precisely on the Rabi 
O ■ frequency, an effect recently called "driven decoherence" . By direct numerical solution of the time- 

dependent Schrodinger equation of the associated many-body system, we scrutinize the different 
mechanisms that may lead to this type of decoherence. Assuming the effects of dissipation to be 
_ -. . negligible (Ti = oo), it is shown that a system of dipolar-coupled spins with - even weak- random 

inhomogeneities is sufficient to explain the salient features of the experimental observations. Some 
experimental examples are given to illustrate the potential of the numerical simulation approach. 
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Decoherence generally occurs when the phase angle associated with a periodic motion is lost due to some interaction 
with exterior noise. In classical mechanics it may apply to classical waves such as sound waves, seismic waves, sea 
waves, whereas in quantum mechanics it applies to the phase angles between the different components of a system 
in quantum superposition. The loss of phase of a quantum system may bring it to its classical regime, raising the 
question of whether and how the classical world may emerge from quantum mechanics. Together with the claim that 
decoherence is also relevant to a variety of other questions ranging from the measurement problem to the arrow of time, 
■ this underlines the important role of decoherence in the foundations of quantum mechanics. It is for all these reasons 
that the analysis of decoherence in quantum systems must make allowance and in particular must distinguish between 
I . decoherence induced by the imperfections of real systems and intrinsic decoherence induced by identified or hidden 
' couplings to the environment. The different sources of decoherence can be classified in two main categories the 
7-H . one-qubit decoherence coming from the coupling of individual qubits with the environment 043 and the multi-qubit 

or pairwise decoherence coming from multiple interactions between pairs of qubits [H-IH 
• i-H , In this paper we take the example of paramagnetic spins because of the quality of the systems which can be 
elaborated (single-crystals) and the possibility, offered by magnetism, to start calculations from first principles. Here, 
the one-qubit decoherence is, in general, associated with phonons and hyperfine couplings P4l2| which are intrinsic 
effects, but also with non-intrinsic effects resulting from weak disorder always present in real systems of finite size: 
inhomogeneous fields, g-factor distributions, and positional distributions. Multiple-qubit decoherence is generally 
due to pairwise dipolar interactions with distant electronic or nuclear qubits, which is an intrinsic mechanism |13| . 
Below, we shall see that, more generally, when pairwise decoherence takes place in the rotating frame, extrinsic 
decoherence becomes crucial by itself and also by amplifying intrinsic decoherence. In particular, by way of some 
examples, it will be shown that the origin of driven decoherence is of the one-qubit type i.e. with multiple possible 
origins (depending on the nature disorder). Even if dominant sources of decoherence may sometimes be identified, the 
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complete description of decoherence and in particular, the discrimination of intrinsic and extrinsic decoherence are 
generally not accessible to experimentalists. This is a major obstacle for the reduction of decoherence, and it holds 
beyond magnetism. We believe that the present, pragmatic approach, should be of great help in common situations 
where intrinsic and extrinsic decoherence mechanisms interoperate. 

Assuming that each type of decoherence has its own "signature" on the Rabi oscillations, we have started a 
systematic study in which the Rabi oscillations of an ensemble of spins are simulated by direct numerical solution of the 
time-dependent Schrodinger equation (TDSE) of the associated many-body system. These simulations are per formed 
using a parallel algorithm implementation based on a massively parallel quantum computer simulator [14]. The 
various mechanisms that may lead to decoherence of Rabi oscillations are successively implemented in Hamiltonians, 
leading to different types of damping, oscillation shapes, non-zero oscillation averages and their evolutions with 
exterior parameters such as the microwave power and the applied static field. The comparison with measured Rabi 
oscillations allows us to scrutinize the different decoherence mechanisms and to understand more basic aspects of 
decoherence, thereby opening a route to search for the optimal - intrinsic and extrinsic - ways to improve coherence 
of Rabi oscillations, i.e. the number of oscillations which is important for all applications. 

The present study is limited to the decoherence of Rabi oscillations, that is the decoherence measured immediately 
after the application of a long microwave pulse. Following an earlier suggestion fl5l - fi"7j |. it was shown that the 
microwave pulse inducing Rabi oscillations is itself an important source of decoherence in all the investigated systems 
("driven decoherence" |18l - l21| ). except when the microwave power is very small, in which case the Rabi frequency is 
also very small. As a consequence the number of Rabi oscillations remains nearly constant, that is one cannot increase 
it by increasing the microwave power. 

This observation can be quantified by comparing the damping time of Rabi oscillations (Rabi decay time Tr) with 
the usual spin-spin relaxation time T%- The theoretical results given in this paper are all exact. Depending on the 
Hamiltonian parameters, the results were obtained analytically (in simple cases) and numerically (in more general 
cases, including dipolar interactions) and covered the large range of possibilities, namely from Tr <C T2 upto Tr w TTi 
when dipolar interactions dominate (in the absence of disorder). 

The systems used to compare the simulations results with experimental data are insulating single-crystals of 
CaW04:Er 3+ , MgO:Mn 2+ , and BDPA (a — 7-bisdiphenylenc-/3-phenylally), a free radical system often used in Elec- 
tron Paramagnetic Resonance (EPR) calibration. The latter is not a diluted system, contrary to the two others, but 
an antiferromagnetic single crystal (identical environments) with a Neel temperature much smaller than the tempera- 
ture at which our measurements are made (between 4K and 300K). These systems have been chosen in particular for 
the differences in their homogeneous/inhomogeneous EPR linewidths. Furthermore, in these systems the relaxation 
time T\ is much larger than T2, as this is often the case in solid state systems. For instance, our experiments yield a 
T\ which is 10 and 40 times larger than the T2 for MgO:Mn 2+ and CaW04:Er 3+ , respectively. Therefore, as a first 
step in the theoretical modeling of these experiments, it is reasonable to neglect the effect of dissipation and focus on 
the decoherence only. 

Rabi oscillations measurements have been performed in a Brukcr Elexsys 680 pulse EPR spectrometer working at 
about / = 9.6 GHz (X-band). Depending on the sample, measurements have been done at room temperature down 
to liquid helium temperature (4K). The static magnetic field has always been chosen to correspond to the middle 
of the EPR line. The experimental procedure is illustrated in Fig. [T] A microwave pulse Pr starts at t = and 
coherently drives the magnetization. At the end of the pulse (r — tn) the magnetization is recorded. Because of 
the dead time (T<jt) of the spectrometer (about 80ns), it is impossible to directly measure the magnetization right 
after the pulse Pr. In this paper, we used two methods for the detection. The first one is simply to record the free 
induction decay (FID) emitted by the system when the microwave field is shut down. This method gives the value 
of the magnetization component M v (t) at the end of the pulse if we take into account two important conditions: i) 
Pr is a non selective pulse (all spins of the line are excited) and under this condition, the FID signal is the Fourier 
transform of the EPR line, ii) The EPR linewidth must be sharp enough. Since the FID is the Fourier transform of 
the EPR spectrum, a linewidth > AG will lead to a decay time of the FID less than 80 ns and the FID will be hidden 
by the dead time of the spectrometer. The second method is used when the EPR line is too broad or if one wants to 
probe the longitudinal magnetization M z (t). In this case another probe sequence has to be used. After the Pr pulse, 
one waits a time T much longer than T2 but smaller than T\ in order such that M v (t) vanishes. After the waiting 
time T, a standard Hahn echo sequence (ir/2 — t — ir — r— echo) is used to measure the longitudinal magnetization. 
In the present paper we do not study (a) the effects of the spin-echo pulses on the measurements and (b) the effect 
of temperature. For (a), this implies that the comparison with theory is through the measured so-called free-decay 
time T 2 * (different from the usual T2) in which a component of the total magnetization is directly measured through 
an induction coil, and for (b) that the measurements of T£ are done at a sufficiently low temperature which is quite 
easy to realize since the T 2 * of BDPA is nearly independent of temperature, and more generally the Rabi time Tr, 
most important in the context of this paper, also. 

The paper is organized as follows. In Section [Til t ne quantum spin model is specified in detail and the simulation 
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FIG. 1: (Color online) Pulse sequence used for Rabi oscillation measurements, a) At t = 0, a microwave pulse Pr coherently 
drives the magnetization. At the end of the pulse (r = rjv) the magnetization component M v (t) is recorded, b) After the Pa 
pulse, one waits a time much longer than T2 but smaller than T\ such that M v (t) has vanished. Then, after the waiting time 
T, a standard Hahn echo sequence of duration rp is used to measure the magnetization component M z (t). The dead time of 
the spectrometer is denoted by r dt . 



procedure is briefly discussed. Our results are presented in Section Hill As there are many different cases to consider, 
to structure the presentation, the results have been grouped according to the kind of randomness, describing for each 
kind (i) the non-interacting case, (ii) the interacting case and (iii) a comparison with experiments if this is possible. 
In Section UVl we present a model of "averaged local Bloch equations", giving a complete, exact description of one- 
qubit decoherence and incorporates multi-qubit decoherence phenomenologically. A summary and outlook is given in 
Section El 



II. MODEL 



We consider a system of L dipolar-coupled spins subject to a static magnetic field along the z-axis and a circular 
polarized microwave perpendicular to the z-axis. The Hamiltonian reads 



fJ-ofJ-B \ " Sj ■ g 3 ■ g fc ■ S fc (Sj ■ gj ■ r 3 fc)(S fc ■ g fc ■ T jk ) 



4tt 



j<k 



(1) 



where Bj(i) = (B'j coscjt, — B'j sinwi, Bq) denotes the external magnetic field, composed of a large static field Bq 
along the z-axis and a circular time-dependent microwave field B'^ (maxj \B'j \ <C minj \Bq\) which may depend on the 
position of the jth spin, represented by the spin-1/2 operators Sj = (SJ , , S?) with eigenvalues ±1/2. The vector 
Tjk connects the positions of spins j and fc. It is assumed that the g-tensor can be written as gj = g e (TL + Agj) where 
the perturbation Agj is a random matrix. 

As usual in the theory of NMR/ESR, we separate the fast rotational motion induced by the large static field Bq 
from the remaining slow motion by a transformation to the reference frame that rotates with an angular frequency 
determined by Bq. Taking the ideal, non- interacting system without fluctuations in the g-tensors as the reference 
system, we define ujq — g e fiBBo and assume from now on that this ideal system is at resonance, that is the microwave 
frequency is given by uj = luq. 

The transformation to the reference frame rotating with angular frequency luq is defined by 



Xrf — 



exp 



ituii 



X exp 



L 



(2) 



where X denotes any combination of spin operators. Transforming Eq. (fT|) to the rotating frame, we find that H^f 



contains contributions that (i) do not depend on time, (ii) have factors 



O ituo 



or e 



-itujQ 



or (iii) have factors e 2j * w ° 
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or e _2tt "°. Contributions that depend on time oscillate very fast (because ojq is large) and, according to standard 
NMR/ESR theory, may be neglected, which we have confirmed for a few cases. The remaining time-independent, 
secular terms yield the Hamiltonian 

H RF = -UJ J2 9 \ 9C S- - fiB9e J2 B 'j 9l 0„ 9 3 S 3 

3=1 9e 1 = 1 9e 

. ^We ^ ^9 z k [l-^ k /r%} 9!g x k ll-3x%/r%}+g y 3 9 y k ll-3y%/r%} 

47F 3<k 9 e T jk A 9e r jk 

From Eq. (|3j), it is clear that variations in have the same effect as local variations (inhomogeneities) in the static 
magnetic field. Inhomogeneities in the microwave field and the variations in (<?J,<7j) are cumulative. Although the 
variations in , g v -, and g* also affect the dipolar interactions, these effects may be difficult to distinguish from the 
effect of the positional disorder of the spins in the solid, in particular if the spins are distributed randomly. Note that 
the total magnetization M z — X^=i &j commutes with the dipolar terms in Eq. Therefore, neglecting the terms 
that oscillate with u>q or 2u>q, the longitudinal magnetization M z (t) = &j ^ s a constant of motion in the absence 

of a microwave field (B'j = 0). 

If all the g's are the same and equal to g e the Hamiltonian Eq. ([3]) reduces to the familiar expression 



h rf = - jr hj s° + £ [1 fast ~ \ (sfss + syst) , . i , 

3=1 j<k 3 k L 



2 

of the Hamiltonian of dipolar-coupled spins in the reference frame that rotates at the resonance frequency ujq 



A. Simulation model 



We now specify the model as it will be used in our computer simulations. We rewrite the Hamiltonian Eq. ([3]) as 



H RF /h = -2ttF £JS| - 2nh p F R £ ( 1 + 0)(2 + g,- ±jj} S * 

3=1 3 = 1 



j< k 3 k 

where we take Fq = ujo/2irh — 9.7 GHz for the Larmor frequency induced by the large static field, Fr = 55.96 MHz 
denotes the Rabi frequency for an isolated spin in a microwave field of 1 mT, we introduce h p as a parameter to control 
the amplitude of the microwave pulse, Do = 51.88 GHz, and we express all distances in A. With this choice of units, it 
is convenient to express frequencies in MHz and time in fj,s. The new dimensionless variables for a — x,y, z and Q 
are defined by g" = g e (l + £") and nsgeBj/h — 27t/i p Fr(1 + Q), respectively. For concreteness, we assume that the 
spins are located on the Si diamond lattice with lattice parameter 5. 43 A (to a first approximation, the choice of the 
lattice is not expected to affect the results). Not every lattice site is occupied by a spin: We denote the concentration 
of spins (number of spins/A 3 ) by n. In experiment n w I0 -4 , . . . , 1 0~ 6 . 

Guided by experimental results, we assume that the distribution of is Lorenztian and independent of a, cut-off 
at £o> and has a width T: 

The reason for introducing the cut-off £0 > is that because the Lorentzian distribution has a very long tail, in practice, 
we may generate £'s such that the corresponding value of g is negative, which may not be physical. Therefore, we use 



= Ttan[(2r - 1) arctan(£ /T)], 



(7) 
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to generate the random variables with distribution Eq. © from uniformly distributed random numbers < r < 1 . 
Likewise, the microwave amplitudes are given by Bj — h p (l + Cj) where the £'s are random numbers with distribution 

"fe»° a,ct m K./7)c;+y eK °- |Cil) ' (8) 

and /ip is the average amplitude of the microwave field. Appendix [3] gives a summary of the model parameters that 
we use in our simulations. 



B. Simulation procedure 



The physical properties of interest, in particular the decay rate cr, = 1/Jr of the Rabi oscillations and the intrinsic 
decay rate C2 = l/?2, can be extracted from the time-dependence of the longitudinal and transverse magnetization, 
respectively, and are defined by 

L L 

<M*(t)> = £<SJ> = £(*(i)|S||*(i)) (9) 
3=1 i=i 

respectively. We compute the time-dependent wave function by solving the TDSE 

ihl\#(t)) = H RF m)) ) (11) 

with Hrf given by Eq. ((5]). Numerically, we solve the TDSE using an unconditionally stable product-formula algo- 
rithm [22] ■ For the largest spin systems, we perform the simulations using a parallel implementation of this algorithm, 
based on a massively parallel quantum computer simulator [l4j ]. Our numerical method strictly conserves the norm 
of the wave function and conserves the energy to any desired precision (limited by the machine precision). 

In analogy with the experimental procedure, we carry out two types of simulations yielding the longitudinal (trans- 
verse) magnetization {M z (t)) {{M x {t))). From Eq. ©, it follows directly that d{M z (t))/dt = -(M«(t)) s hence 
(M z (t)) is directly related to (M y (t)) measured in experiment. We prepare the spin system, that is the state |^(0)}, 
such that all spins are aligned along the z (x)-axis. Then, for a fixed value of the microwave amplitude h p (h p = 0) 
and a particular realization of the random variables £J, £j, £J, Q and the distribution of the spins on the lattice, we 
solve the TDSE and compute Eqs. (|9"1)- (fTU|) . This procedure is then repeated several times with different realizations 
of random variables and distributions of spins. Finally we compute the average of Eqs. (Ii?])- (I10I) over all these realiza- 
tions and analyze its time dependence by fitting a simple, damped sinusoidal function to the simulation data. This 
then yields the decay rate cr = 1/Tr (intrinsic decoherence rate C2 = I/I2) °f the Rabi oscillations. 



III. RESULTS 



In the subsections that follow, we consider the various sources of decoherence separately. We also study the interplay 
of intrinsic decoherence due to e.g. pairwise interactions and extrinsic decoherence due to e.g. single spins driven by 
external magnetic fields when different spins have different environments (different couplings to static and microwave 
fields). The averaging over different spins leads to decoherence, that is phase destruction of the electromagnetic waves 
generated by the spins. These two types of decoherence lead to the observed damping of Rabi oscillations which 
takes place through energy exchange between the spin system and the applied microwave field. Below we show that 
energy dissipation from the spin-bath to the electromagnetic bath is sufficient to explain the experimental results on 
the Rabi decay time. This is the reason why we neglect, in the present paper, the dissipation effect of phonons, (our 
spin-lattice relaxation time T\ is infinite). Note that if we turn off the microwave field, the longitudinal component of 
the magnetization commutes with the Hamiltonian Eq. ([5]) and hence does not change with time at all, showing that 
energy exchange with the electromagnetic bath is essential. In the following sections we give two examples in which 
energy flows from the electromagnetic bath to the spin-bath and from the spin-bath to the electromagnetic bath. 
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FIG. 2: (Color online) The Rabi oscillations of the longitudinal magnetization as obtained by solving the TDSE for the 
Hamiltonian Eq. (T5J for 26 dipolar-coupled spins for two different concentrations n, without random fluctuations in the g- 
factors (r = 0) and on the microwave amplitude (7 = 0). Top left: n = 10 -3 ; Top right and bottom: n = 10~ 4 . The solid lines 
represents the envelope (a ± be~ CR *)/2 of the function (a + be~ CRt cos 2wft) /2 that was fitted to the data. 

A. Fixed (/-factors and homogeneous fields 

In the absence of randomness on the g-factors or on the microwave amplitude, the Hamiltonian is given by Eq. ([5]) 
with £ x = £ y = = C = 0. 



1. Non-interacting spins 



For non-interacting spins, we can drop the spin label and write the Hamiltonian (in the rotating frame) for a single 
spin as 



H RF /h = -2nh p F R S x . 
The time evolution of the longitudinal magnetization takes the simple form 



(12) 



(¥(f)|S'|*(t)> = -cosft R i, 



(13) 



showing that the ^-component of the spin performs undamped Rabi oscillations with angular frequency Or = 2Trh p F^i. 
Therefore Tr = 00. Furthermore, the transverse magnetization is conserved and therefore T2 = 00. Summarizing, in 
the absence of randomness and dipole-dipole interactions, we have 



T R = 00 , T 2 — 00. 



(14) 
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t (us) t (us) 

FIG. 3: (Color online) Time evolution of the transverse magnetization as obtained by solving the TDSE for the Hamiltonian 
Eq. ([5j) for 26 dipolar-coupled spins without random fluctuations in the (/-factors (r = 0) and on the microwave amplitude 
(7 = 0). Left: n = 10~ 3 ; Right: n = 10 -4 . The non-zero value of a is due to the statistical noise and the relatively short time 
interval used for the fit. The solid line represents the function (a + be~ C2t )/2 that was fitted to the data. 

2. Dipolar- coupled spins 

In Figs. [2] and [3l we present simulation results for the longitudinal and transverse magnetization, respectively, as 
obtained by averaging the solutions of the TDSE over ten different distributions of 26 dipolar-coupled spins on the 
lattice. Our simulation results, many of them not shown, lead us to the following conclusions: 

• For both concentrations n = 1CP 3 and n — 10 -4 and for microwave amplitudes h p — 0.5,1,2, the Rabi 
oscillations decay exponentially. Indeed, the fits are good, as indicated by the small differences between the 
Rabi frequency (Fq = 55.96) and the values of / obtained by the fitting procedure. 

• The decay rate cr, = 1/Tr increases with n, with a slope of approximately 1.7 (data not shown). 

• Within the statistical fluctuations resulting from the random distribution of the spins on the lattice, cr = 1/Tr 
does not depend on the microwave amplitude h p but strongly depends on the concentration n. 

• Simulations (data not shown) for n = 0.25 x 10 -4 , . . . , 10~ 3 indicate that T 2 <x n, as expected theoretically. 

• The simulation data suggest that c 2 = I/T2 > cr = 1/Tr. 

Summarizing, in the absence of local randomness but in the presence of dipole-dipole interactions, we have 

T R = T R {n)>T 2 =T 2 (n). (15) 

3. Experimental results: BDPA 

We now compare these theoretical predictions to experiments performed on a single crystal of BDPA (a — 7- 
bisdiphenylene-/3-phenylally). With a linewidth of 0.09 mT, this system is quite homogeneous with a very narrow 
distribution of the g-factors. Moreover, the sample used was very tiny such that we may consider the microwave to 
be homogeneous inside the sample. 

Results are presented in Fig. 2J They show an example of Rabi oscillations obtained from FID measurements. The 
Rabi oscillations fit very well to 

M v (t) = A a sm(n R t + <j))cxp{-t/T R ) + M v {oo), (16) 

for all microwave powers. The obtained Rabi decay time Tr is clearly independent of the amplitude of the microwave 
field, as predicted by the model when T = 7 = 0. It is also very close to T 2 *, the FID decay time given by the Fourier 
transform of the EPR linewidth. This is also in agreement with predictions when T = 7 = and Dq ^ 0, T 2 * being a 
coherence time fully equivalent to T 2 . The discrepancy between Tr (~ 140 ns) and T 2 (= 128 ns) is due to a small 
inhomogeneous broadening (about 10%). 
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FIG. 4: (Color online) a) Decay time Tr of the Rabi oscillation of BDPA as a function of the amplitude of microwave field. 
Dashed line: Value of TJ the decay time of the FID. b) Example of the FID signal. Tdt- Dead time of the spectrometer. Black 
line (red line): In-phase (out-of-phase) FID signal recorded by the detector. Blue line: Best fit to the exponential decay, c) 
Example of Rabi oscillations. The red line is the best fit to M y (t) = Ao sm(ilnt + <j>) exp(— t/Tu) + M v (oo), giving the Rabi 
decay time Tr and a non-zero offset. The nonzero value of M y (oo) is due to dissipation effects [23], collectively described by 
the relaxation time Xi, which are not included in the microscopic model considered in the present paper. Measurements were 
carried out at room temperature. 



B. Randomness in the microwave amplitude only 



In the case of randomness in the microwave amplitude only, the Hamiltonian is given by Eq. ([5]) with £ x — t; y — 
£ z = 0. Such a randomness is inherent to finite size cavities and becomes smaller as the size of the sample relative to 
the size of the cavity is reduced. 



1. Non-interacting spins 



For non-interacting spins (Do = 0), we can readily compute the average over the distribution of Q analytically if 
we neglect the cut-off of the Lorentzian distribution. As all spins are equivalent, we may drop the spin index j and 
we obtain 

showing that the Rabi oscillations decay exponentially and that the decay time of the Rabi oscillations is given by 
1/Tr = 7^r. Furthermore, the transverse magnetization is conserved and therefore T2 — 00. Summarizing, in the 
presence of randomness in the microwave field only and in the absence of dipole-dipole interactions, we have 

1/T R = 7 Q R > 1/T 2 = 0, (18) 

showing that the decay rate of the Rabi oscillations increases linearly with the microwave amplitude h p whereas T2 
remains infinite. This is easy to understand: T2 is infinite due to the lack of pairwise intrinsic decoherence whereas 
destructive interference associated with weak positional randomness in h p (the microwave field) leads to a reduction 
of Tr when h p increases (one-qubit decoherence). 
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FIG. 5: (Color online) Simulation results as obtained by solving the TDSE for the Hamiltonian Eq. (JS| for 12 spins (concentration 
n — 10 -4 ) that interact via dipole-dipole interaction, without random fluctuations in the p-factors (r = 0) but with random 
fluctuations on the microwave amplitude (7 = 0.01). The results represent the average of 100 different realizations of 12-spin 
systems. Top left to middle right: Longitudinal magnetization for different values of h p . The solid line represents the envelope 
(a ± be~ CR )/2 of the function (a + fee _CRt cos 2n ft) /2 that was fitted to the data. Bottom left: Transverse magnetization in 
the absence of the microwave field (h p = 0). The solid line represents the function (a + be~ C2t )/2 that was fitted to the data. 
Bottom right: Bullets show the inverse relaxation time cr = 1/Tr as a function of the microwave amplitude h p . The dashed 
line connecting the bullets is a guide to the eye only. A linear fit to the simulation data yields cr = 1/Tr ss 3.69h p + 1.82 
and is shown by the solid line. The horizontal line represents the value of C2 = I/2T2 as 1.95, estimated from the data of the 
transverse magnetization (see bottom left). 
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FIG. 6: (Color online) Decay time Tr of Rabi oscillations in CaWCU as a function of the microwave field amplitude (corre- 
sponding to the Rabi frequency Ojj) for two concentrations of the Er 3+ spins. The static field Ho is parallel to the c-axes of 
the cristal and the temperature T~4 K. Dashed lines: Fit to 1/Tr = 1/Tro + f3£lu/2n Measurements were carried out at 4K. 




FIG. 7: (Color online) Rabi oscillations of MgO:Mn 2+ (0.001%) for two sample sizes. Measurements were carried out at room 
temperature. 



2. Interacting spins: dipole-dipole interaction 

In Fig. [SI we present simulation results for systems of 12 spins with dipole-dipole interaction and randomness in 
hp, as obtained by averaging over 100 different realizations, meaning 100 different distributions of the 12 spins on the 
lattice. The four upper panels of Fig. [5] show results for the longitudinal magnetization (M z (t)). 

Rabi oscillations are damped but have zero offset. The inverse Rabi time 1/Tr = cr, deduced from sinusoidal fits, 
increases linearly with the microwave field, that is with the Rabi frequency Qr (bottom right). Its value at h p = is 
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to good accuracy equal to I/2T2 (ao ~ 1.82 for n — 10~ 4 ). The slope a\ ~ 3.69/Fr is related to the matrix of the 
gyromagnetic factor and to the root mean square of local fields resulting from the randomness in the microwave field. 
These results, specific to a h p distribution, agree qualitatively with recently published results of the damping time of 
Rabi oscillations in the limit of a large inhomogenous linewith [24{ . 

The results for the transverse magnetization (M x (t)) in the absence of microwaves (h p = 0) are presented in the 
bottom left panel of Fig. [5j It clearly decays exponentially, as this is the case with the longitudinal magnetization. 
Summarizing, from Fig. [5] we conclude that in the presence of randomness in the microwave field and of dipole-dipolc 
interactions, we have 

cr = 1/Tr « aifiR, + ao , a w 1/2T 2 . (19) 

Here, pairwise decoherence affects Ti which is now finite (and proportional to 1/n as in the case without randomness, 
see Section [III Ap and randomness in microwave amplitude h p affects Tr which is essentially proportional to l/h p at 
large h p . As Tr < T2, we can say that, in this case, energy flows from the spin-bath to the electromagnetic bath, 
leading to energy dissipation in the spin-bath. 

3. Experimental results: CaWO i :Ef i+ and MgO:Mn 2+ . 

In order to show the effect of concentration on Rabi damping we measure two samples of CaWC>4:Er 3+ with 
Erbium concentration 0.01% and 0.001%, respectively. The two samples have nearly the same shape, keeping the 
inhomogeneity of microwave field constant. To remove the effects of zero microwave field decay (that is T2 due to 
multi-spin or pairwise decoherence) we plot 1/Tr — 1/Tro where Trq ~ I/2T2 is the decay time at zero microwave 
field. The results are presented in Fig. The inverse Rabi decay time fits very well to 1/Tr = I /Trq + /3£lft/27r, 
where /3 is a fitting parameter. From Fig. [6l it is clear that the Rabi-decay time Tr decreases with the concentration 
n, in concert with the simulation results. 

Evidence of the effect of microwave field inhomogeneity on the Rabi oscillation decay has been recently given for 
a sample of Cr:CaW04 [251 ]. To provide further evidence, we took a sample of MgO doped with about 0.001% with 
Mn 2+ and cut the sample into a large (3.6 x 5 x 0.5 mm 3 ) and small (1 x 1 x 0.5 mm 3 ) piece. At this extremely low 
concentration, the dipole-dipole interaction effect on the Rabi decay is negligible, hence disorder essentially due to 
the microwave field inhomogeneity inside these samples will be different. Fig. [7] shows the Rabi oscillations for these 
two samples. All parameters (microwave power, temperature, crystal orientation) are the same for the measurements 
on these two samples. The effect of the inhomogeneity of the microwave field on the Rabi decay time is clearly seen 
as the damping in the large sample (red line) is almost two times larger than the one in the small sample (black line) . 

C. Randomness in the g- factors only 

We assume that there are no random fluctuations in the amplitude of the microwave pulse and that the g-factors 
fluctuate randomly from spin to spin. This effect is generally due to weak crystal distortions, imperfections, leading 
to small variations of crystal-field parameters. 

1. Randomness in g z : Non-interacting spins 

In this case, the Hamiltonian is given by Eq. ([5} with Q = £ x = £ v = Dq = 0. As we then have a system of 
independent spins, we may drop the spin index j. In the case that initially, all the spins are aligned along the z-axis, 
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In the case that initially, all the spins are aligned along the x-axis, we find 
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Recall that we calculate the transverse magnetization for the case that initially, all spins are aligned along the x-axis. 
In order to obtain the expressions in terms of elementary functions, we have ignored the cut-off of the Lorentzian 
distribution. We can check that for T = 0, Eq. (|2"0"|) and Eq. (f2"Tj) reduce to 
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{S*(t)) = -cosQ R t 



(S"(t)) = \, 



while for h p = 0, we find 



(S*(t)) = - 



(S*(t)) = - e - 2 " trF <> 



(22) 



(23) 



in agreement with the expressions that can be derived directly, without any averaging procedure. From Eq. (|23p . it 
follows that I/T2 = 2ttTFq. For finite T, Rabi oscillations are present only if h p Fn > TFq in both longitudinal and 
transverse cases. 

In Fig. [Sf left), we present a typical result for the time dependence of the longitudinal magnetization with <? 2 -factor 
distribution (only) suggesting that the time-averaged longitudinal magnetization is non-zero, in concert with the 
analytical expressions 



lim i f (S*(t))dt = 0, 
lim - fwW)dt = 0, 



lim - , 

T->oo T J 



(S z (t))dt = 



1 



TF 



2 h p F R + T F 



(24) 
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FIG. 8: (Color online) Left: Time evolution of the longitudinal magnetization as obtained by numerical evaluation of Eq. (|20p 
(or the solution of the TDSE for Do = £ x = f v — C = f° r F = 0.001), that is for the case that there are random fluctuations in 
the <? z -factor only. The solid line represents the envelope (a ± be~ CRt )/2 of the function (a + be~ CRt cos27r ft) /2 that was fitted 
to the data. Right: The decay rate cr as a function of h p , obtained by fitting (a + be~ CRt cos2irft)/2 to the time-dependent 
data. The solid line shows that the function a' \n(h p ) + b' with a' = —0.35 and b' — 1.32 fits the data reasonably well. The 
dashed line connecting the bullets is a guide to the eye only. 



The reason for this positive offset is simple: Any non-zero field in the z-direction tilts the plane of the Rabi oscillations 
away from the (y, z)-plane, introducing a small precession about the tilted axis superimposed on the Rabi nutation, 
leading to a positive long-time average. This non-zero offset effect is significant because, as we will see later, it is a 
unique signature of the presence of random fluctuations in the g z factor or, cquivalcntly, of the inhomogeneity of the 
static magnetic field. We emphasize that this non-zero offset is due to randomness and not due to dissipation, as the 
present paper considers the case of T\ — oo only. 

Similarly, in the case that all spins are initially along the x-direction, the long-time average of the transverse 
magnetization is given by 

jsJ f ™ - \l^kTrK- < 25) 

the long-time averages of the two other components being zero. Unlike in the case of the longitudinal magnetization, 
in the regime where the transverse magnetization shows oscillations {h p F^ > TFq), the transverse magnetization 
reaches its asymptotic value Eq. f|25[) already after a few oscillations (data not shown). 

From Eq. (f2"U|) , it is clear that we cannot expect the amplitude of the Rabi oscillations to decay exponentially in 
a strict sense. Nevertheless, the data fits well to a function of the form (a + 6e~ CR * cos27r/t)/2. The decay rate cr, 
shown in Fig. right), decreases with increasing microwave amplitude h p . It seems to diverge when h p — >• but this 
is never observed in experiment. 

This decrease is a second characteristic feature of the presence of random fluctuations in the g z factor or, equivalently, 
of the inhomogeneity of the static magnetic field. 



2. Randomness in g x and g y : Non-interacting spins 
In this case, the Hamiltonian is given by Eq. with Q = £ z = Dq = and we have 

WW = \ J °°co S [o a t(i+(r+e)/2)]p(rMarfr^- (26) 

Taking the cut-off £o to be infinity we obtain 

i r+°° r i 

(S'(t)) = ^ J cos [fi R t(l + 0] dC = - e -°* rt cosftRt. (27) 

Thus, we conclude that if there is randomness in g x and g v only, the Rabi oscillations will decay exponentially with a 
rate proportional to Qr = 2ixh v Fji- In the absence of the microwave field, the transverse magnetization is a constant 
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FIG. 9: (Color online) Top left: Time evolution of the longitudinal magnetization as obtained by numerical solution of the 
TDSE (see Eq. pip ) for the case that there are random fluctuations in all three gr-factors only (r = T x — F y = F z — 0.001 
and £ = -Do = 0). The solid line represents the envelope (a ± be~ CRt )/2 of the function (a + be~ CRt cos 2irft)/2 that was fitted 
to the data. Top right: Same as top righ, except that h p = 10 instead of h p — 0.5. Bottom left: Transverse magnetization 
in the absence of the microwave field (h p = 0). The decay rate C2 = 60.19 is in excellent agreement with the analytical result 
cr, = 27iTFo = 60.95 predicted by Eq. (J23J) . The solid line represents the function (a + be~ C2t )/2 that was fitted to the data. 
Bottom right: The inverse relaxation time cr as a function of the microwave amplitude h p for F x — F y = F x — 0.001 and 
F z = 0.001 (bullets), F z = 0.002 (squares), F z = 0.003 (triangles). The solid line represents the linear fit to the F z = 0.001 
data. The dashed lines are guides to the eye only. The number of spins in these calculations is 10000. 

of motion and hence T2 = 00. Summarizing, in the presence of randomness in g x and g y only and in the absence of 
dipole-dipole interactions, we have 

1/T R = m R > 1/T 2 = 0, (28) 

showing that the decay rate of the Rabi oscillations increases linearly with the microwave amplitude h p . In fact, 
Eq. (|28|) is the same as Eq. (fT5|) with 7 replaced by T. Thus, we conclude that randomness in g x and g y has the 
same effect as randomness in the amplitude of the microwave field: The Rabi oscillations decay exponentially, with 
a decay rate that increases linearly with Qr = 2Tth p FR. In both cases, decoherence results from a loss of phase of 
superposed radiation emitted by spins in nutation leading, as a consequence, to energy transfer from the spin-bath 
to the electromagnetic bath. Clearly enough such dissipation does not involve the usual relaxation time T\ due to 
dissipation by phonons. This case is very different from the one of e.g. superconducting qubits where decoherence is 
dominated by T\ process, as shown for example in Ref. (26|. 

3. Randomness in g x , g y and g z : Non-interacting spins 

In this case, the Hamiltonian is given by Eq. §5§ with Q = Dq = 0. In Fig. E^top), we present a typical result for 
the time dependence of the longitudinal magnetization. It is seen that the time-averaged longitudinal magnetization 
is non-zero, signaling the presence of fluctuations in g z (see Section UlI C 1[) . Also clearly visible is the increase of the 
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FIG. 10: (Color online) Decay time of the Rabi oscillations of MgO:Mn 2+ (0.001%) as a function of the microwave field 
amplitude (Rabi frequency Qr) for two samples of different sizes. Measurements were carried out at room temperature. 

decay rate cr of the Rabi oscillations with increasing microwave amplitude h p , a signal of the presence of fluctuations 
in (g x ,gy) (see Section [III C 2p . Note that there is no obvious relation between the decay rate of the transverse 
magnetization (c2 ~ 60, see Fig.^bottom left)) and the values of the decay rate cr at the smallest values of h p shown 
in Fig. [^bottom right). 

From the results of Sections MI C II and MI C 21 we may expect that the decay rate cr shows a crossover from the 
regime in which the fluctuations on g z dominate (cr decreases with increasing h p ) and a regime in which the fluctu- 
ations on (g x ,gy) dominate (cr increases linearly with h p ). This is borne out by the data presented in Fig. [^bottom 
right) where we show the combined effect of the two different sources of decoherence, the widths of the Lorenztian 
distributions for the longitudinal (g z , T z ) and transverse ( (g x ,gy), F x = T y ) fluctuations being varied independently. 

4- Experimental results: MgO:Mn 2+ 

The combined effect of a distribution in the g-factors and inhomogeneities in the microwave amplitude are shown 
in experiments performed on single crystalline films of MgO:Mn 2+ , see Fig. 1101 where the measured Rabi dacay time 
is plotted versus the Rabi frequency. The Mn 2+ dilution is such that dipolar interactions are negligible. Due to 
weak but sizable distributions of Mn 2+ local environments, we expect non-negligible and similar distributions of 
the three g-factor components. For small microwave amplitudes, the distribution in the t^-factor gives the dominant, 
nearly constant contribution to the Rabi decay time, which compares well with Fig. ^bottom right). As the microwave 
amplitude increases, the inhomogeneities associated with transverse components take over and 1/Tr increases linearly 
on the log-log scale. Note that the slope of one-half differs from the slope one that we have for the model considered 
in this paper. This is because of the peculiarity of the experimental system where nutation takes place coherently over 
five equidistant levels of the material, an aspect that will be considered in the future. At present, we are interested in 
showing that the departure from the 1/Tr plateau takes place more rapidly with the larger sample as expected when 
the effect of microwave inhomogeneities dominates over the one of g-factor distributions. 

5. Dipolar- coupled spins 

In Fig. lllf top and middle), we present simulation results for systems of 26 spins with dipole-dipole interaction, (with 
different concentrations n), with random fluctuations in the three g-factors and uniform microwave field amplitude. 
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FIG. 11: (Color online) The Rabi oscillations of the longitudinal magnetization as obtained by solving the TDSE with the Hamil- 
tonian Eq. (|5| for 26 spins that interact via dipole-dipole interaction, for different concentrations n, with random fluctuations in 
the three gr-factors (r = 0.001) and without random fluctuations in the microwave amplitude (7 = 0). Top left: n = 10" 3 ; Top 
right to middle right: n = 10~ 4 . The solid line represents the envelope (a ± be~ CRt )/2 of the function (a + 6e _CR * cos27r/t)/2 
that was fitted to the data. Bottom: Time evolution of the transverse magnetization for n = 10 -3 (left) and n = 10~ 4 (right). 
The solid line represents the function (a + be~° 2t )/2 that was fitted to the data. 



These results are obtained by averaging over ten different realizations, meaning ten different distributions of the 26 
spins on the lattice. The striking signature of the presence of fluctuations in g z , namely the non-zero long-time average 
of the longitudinal magnetization, remains untouched by the effects of the dipolar interactions. For the values of h p 
shown in Fig. ITTT top left to middle right), the dependence of the decay rate cr is essentially the same as if the dipolar 
interactions were absent (see Fig [^bottom right)). For large h p (data not shown), the decay rate cr linearly increases 
with hp. Comparing Fig. ITTT bottom left) with Fig. [TTT bottom right), it follows that the value of the decay rate of 
the transverse magnetization is nearly independent of the concentration, hence cannot be attributed to the presence 
of dipolar interactions but is mainly due to the presence of fluctuations in g z . 
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FIG. 12: (Color online) The Rabi oscillations as obtained from the numerical solution of the TDSE (see Eq. (JTTJ) for Do = 0, 
7 = 0.01 and V = 0.001, that is for the case that there are random fluctuations in both the microwave field and in the (/-factors. 
The solid line represents the envelope (a ± be~ CRt )/2 of the function (a + 6e _CR * cos27r ft) /2 that was fitted to the data. The 
number of spins in these calculations is 10000. 



D. Randomness in the g- factors and the microwave amplitude 

1. Non-interacting spins 

In Fig. [T^J we present a few representative results for the case that there are random fluctuations in both the 
microwave amplitude and in the g-factors, as obtained by solving the TDSE for the Hamiltonian Eq. ((5]) with Dq = 0. 
In essence, the results are very similar to those of the case where there are fluctuations in all three g-factors only. 
This is easy to understand from Eq. ([5]): Fluctuations in (g Xl g y ) or (exclusive) in the microwave amplitude have the 
same effect on the decay of the Rabi oscillations. With both types of fluctuations present, our numerical results show 
that this contribution does not significantly alter the dependence of cr, on h p . 

As before, the presence of fluctuations in g z (see Section MI C II) is signaled by the time-averaged longitudinal 
magnetization being non-zero and by a contribution to the decay rate C2 of the transverse magnetization, which is in 
excellent agreement with the analytical result C2 = 2ttTFq predicted by Eq. (I2U1) (data not shown). Thus, in this case, 
we obviously have c-i > cr which is the same as Tr > T-i where T-i is reduced by the fluctuations in g z . 



2. Dipolar- coupled spins 

In Fig. 1131 we present simulation results for systems of 12 spins with dipole-dipole interaction, as obtained by 
averaging the solution of the TDSE over 100 different distributions of the 12 spins on the lattice, for the case that 
there are random fluctuations in the microwave amplitude and in all three (/-factors. 

The four upper panels of Fig. [T2] show results for the longitudinal magnetization. The decay of the longitudinal 
magnetization is exponential to good approximation. The signature of the presence of fluctuations in g z , namely 
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FIG. 13: (Color online) Simulation results as obtained by solving the TDSE for the Hamiltonian Eq. (J3J) for 12 spins (concen- 
tration n = 10~ 4 ) that interact via dipole-dipole coupling, with random fluctuations in the g-factors (r = 0.001), and with 
random fluctuations in the microwave amplitude (7 = 0.01). Top left to middle right: Longitudinal magnetization showing Rabi 
oscillations. The solid line represents the envelope (a ± 6e _CR ')/2 of the function (a + be~ CR cos 2nft)/2 that was fitted to the 
data. Bottom left: Bullets show the inverse relaxation time cr = 1/Tr as a function of the microwave amplitude h p . The dashed 
line connecting the bullets is a guide to the eye only. A linear fit to the simulation data yields cr = 1/Tr. » 3.98h p + 4.43 
and is shown by the solid line. The horizontal line represents the value of I/2T2 ~ 1.95, estimated from the data of the 
transverse magnetization in the absence of random fluctuations in the g-factors and on the microwave amplitude (see Fig. [5|. 
Bottom right: Transverse magnetization in the absence of the microwave field (h p = 0). The solid line represents the function 
(a + be~ C2t )/2 that was fitted to the data. The decay rate C2 contains contributions from the dipolar interactions and, most 
importantly, from the random fluctuations in g z . 
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the non-zero long-time average of the longitudinal magnetization is clearly visible. For the values of h p shown in 
Fig- fTST bottom left), the linear dependence of the decay rate cr is essentially the same as if the dipolar interactions 
were absent (see Fig ^bottom right)). 

A linear fit to the data of cr yields lim^ _>o cr ~ 4.43. This value should be contrasted with the result C2 f» 63.74 
for the transverse magnetization in the absence of microwaves (h p = 0) (see Fig. [TBT bottom right)). Such a large 
C2 (small T2) resulting from both dipolar interactions and fluctuations on all g- factors is effectively caused by the 
effect of (^-fluctuations, in concert with the results shown in Fig. lllf bottom') that demonstrate that the concentration 
dependence is weak, implying that the effect of the dipolar interactions is small compared to that of the presence of 
fluctuations in g z . 

According to theory, the total decay rate of the transverse magnetization is the sum of the decay rates due to the 
dipolar interactions only and the combined decay rate due to field inhomogeneities only. From Fig. [5l the former is 
given by C2 ~ 3.90. In the absence of dipolar interactions, the latter is given by C2 = 2irTFo — 60.95 (see Section UlI C 3[ 
and Fig. [HI yielding C2 ~ 60.19 for 7 = 0). Therefore, we have 4 otal 64.85, in very good agreement with the value 
c 2 = 63.74 extracted from the simulation (see bottom right panel of Fig. 

IV. PHENOMENOLOGICAL MODEL 

The simulations of the dipolar-coupled spin systems are rather expensive in terms of computational resources. For 
instance, one simulation of a single realization of a 26-spin system takes about 20 hours, using 512 CPUs on an IBM 
BlueGene/P. Such relatively expensive simulations are necessary to disentangle the various mechanisms that may 
cause decoherence but are not useful as a daily tool for analyzing experiments. Therefore, it is of interest to examine 
the possibility whether a simple phenomenological model can capture the essence of the physics of the full microscopic 
model. Based on our results, presented in Section IIII[ we propose to use a single-spin model to which we artificially 
add a dephasing/relaxation mechanism. 

Specifically, we propose that the Heisenberg equation of motion (in the rotating frame) of the expectation values 
of the spin-components is modified according to 

g ( -I/T2 27^**0 \ 

or(S(t)> = -27t£*F -I/T2 Trh p (l + ()(2 + e+e j )F R \ (S(t)), (29) 

where we adopt the same notation as the one used in Section III Al The phenomenological aspect enters in the 
introduction of the decay times T± and T2 . 

Equation (|29|) has the same structure as the Bloch equation but there is a conceptual difference and a practical 
consequence. The former comes from the introduction of g-factor and microwave field amplitude distributions and the 
latter offers the possibility to calculate numerically the effects of one-spin decoherence to a high degree of accuracy. 
As we showed in this paper, one-spin decoherence plays an essential role when several qubits act at the same time. It 
is then natural to start from the well-known equation of motion of a spin S = 1/2, add disorder through distribution 
probabilities (here of g-factors and microwave field amplitude) and average over the solutions. This leads to the exact 
knowledge of corresponding one-spin decoherence, namely to Eq. (f2"9")) without the T\ and T2 terms. If we now want 
to make a link with the Bloch equations we have just to add the phenomenological damping times T± and T2 as it is 
done in the original Bloch equations. The difference between Eq. (|29p and the original Bloch equations is that in the 
latter T\ and T2 include all damping contributions i.e. many-spin and one-spin damping, whereas in the former T\ 
and T2 include many-spins damping only, one-spin damping being calculated exactly. 

Before assessing the usefulness of Eq. (f2"9"| by comparing its results to the numerical solution of the TDSE of the 
interacting spin system, it is instructive to analyze the case £ x = £ v = = £ = T\ = 0. Then the solution of Eq. (|2"9"|) 
reads 

(S x (t)) = e-'/ T2 (S*(0)} 

(S z (t)) = e-'/ 2T2 cos (2irh p F K ^l - (l/4nh p F R T 2 )^j (S z (0)), (30) 

where, for simplicity, we have assumed that (S v (0)) = 0. From Eq. (l30|) it follows that the transverse and longitudinal 
magnetization decays exponentially with a relaxation time T2 and 2T2, respectively. In other words, in the absence 
of randomness and for T% = 0, Eq. (I2TH) predicts a factor of two between the relaxation time of the Rabi oscillations 
and the relaxation time of the transverse magnetization, in qualitative (and almost quantitative) agreement with our 
simulation results of dipolar-coupled spin-1/2 systems with randomness. Thus, model Eq. (|29l) may give a simple 
explanation why in our simulations, we find that extrapolation of cr to h p = gives, in the presence of dipolar 
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FIG. 14: (Color online) The Rabi oscillations of the longitudinal magnetization as obtained from the numerical solution of the 
phenomenological model Eq. (|29[) with Ti = oo, using 10000 realizations of the random variables £ x , £ v , £ z , and Left: V = 0, 
7 = 0.01, 1/Ti = and T 2 = 3.0, compare with Fig. EJtop right); Right: T = 0.001, 7 = 0, 1/Ti = and T 2 = 3.6, compare 
with Fig. [Til middle right). The solid line represents the envelope (a ± 6e~ CR ')/2 of the function (a + be~ CRt cos 2irft)/2 that 
was fitted to the data. 



interactions, precisely C2/2 if there is no distribution of ^-factors (Do 7^ 0, T — 0) and a value larger than C2 if there 
is a distribution of g-factors (r > 0). 

If we put £ x = £, v = £ z = ( = 0, which in principle we should do if we strictly adopt the Bloch-equations approach, 
we can never recover the linear dependence of the decay rate 1/Tr on the microwave amplitude h p . However, if we 
average over the £'s and/or £ and put T2 = 00, the results are the same as those obtained from the direct solution of 
the TDSE of the spin-1/2 system. 

In appendix IBI we give a simple, robust, unconditionally stable algorithm [23] to solve Eq. (|2"5|) . In Fig. RH1 we 
present some representative results. We used the same parameters for T, 7 and h p and changed the phenomenological 
parameter T2 until we found a fair match with the data of the corresponding interacting system. Taking into account 
that we did not attempt to make a best fit to these data, the agreement is excellent. In both cases shown in Fig. [HI 
(and in many others cases not shown), this simple procedure seems to work quite well. This suggests that the simple 
model Eq. (|29p may be very useful for the analysis of experimental data, including the effects of the pulse sequence 
and pulse shapes, effects that are rather expensive to analyze using the large-scale simulation approach adopted in 
the present paper. 



V. SUMMARY AND OUTLOOK 



The main results of this paper may be summarized as follows: 

• The non-interacting spin model can account for the ^-dependence of the decay of the Rabi oscillations if we 
introduce randomness in the g- factors (all three) and/or in the amplitude of the microwave field. In the case 
of g z randomness, the long-time average of the longitudinal magnetization deviates from zero. This deviation 
increases as the Rabi frequency decreases and reaches its maximum (1/2) when hpFft/TFo — > 0. The effect 
of the g z distribution on the value of cr at zero microwave field (h p = 0) is simply related to the value of 
(M z (t = 00)), suggesting that this decoherence effect comes from the combination of different spin precessions 
about the z-axes and the nutational motion of spins. 

• The dipolar-coupled spin system without randomness in all three ^-factors and without randomness in the 
amplitude of the microwave field, cannot account for the VLr- dependence of the Rabi oscillation decay rate, 
observed in experiment. The decay rate of the Rabi oscillations increases as the concentration of magnetic 
moments increases, as one naively would expect. 

• The dipolar-coupled spin system without randomness in g z but with randomness in the amplitude of the mi- 
crowave field and/or randomness in (g x ,gy), can account for the /ip-dependence of the Rabi oscillation decay 
rate and also for the concentration-dependence of this decay rate, just as in the case of non-interacting spins. 

• The dipolar-coupled spin system with randomness in all three the ^-factors and with or without randomness 
in the amplitude of the microwave field, can account for the /ip-dependence of the Rabi oscillation decay rate 
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and also for the concentration-dependence of this decay rate. A salient feature of the presence of fluctuations 
on g z (or, equivalently on inhomogeneities in the static field) is that the long-time average of the longitudinal 
magnetization deviates from zero, as in the case of non-interacting spins. 

For future work, we want to mention that the effects on the decay of the Rabi oscillations of the measurement 
by the spin-echo pulses themselves may be studied by the simple phenomenological model described in Section IIVI 
Among other aspects, not touched upon in the present study, are the case where motional narrowing is important [28j 
or where dipolar interactions are strong enough to induce decoherence by magnons, as recently shown in the Fes single 
molecular magnet [2^ | . These cases can be treated by the simulation approach adopted in this paper and we plan to 
report on the results of such simulations in the near future. 
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Appendix A: Overview of the model parameters 

For convenience, we list the parameters of our model: 

• The Larmor frequency i*o = wo/27rfi, = 9700 [MHz] which is fixed. 

• The Rabi frequency at a microwave amplitude of 1 mT is = 55.96 [MHz] which is fixed. 

• The amplitude of the microwave pulse, controlled by the parameter h p . By convention, if h p = 1, a single 
isolated spin will perform Rabi oscillations with a frequency of Fr = 55.96 [MHz]. The Rabi pulsation in the 
microwave field h p is SIr — 2?t FrH v . 

• The width 7 of the Lorentzian distribution of the random fluctuations of the amplitude of the microwave pulse 
h p . 

• The width T of the Lorentzian distribution of the random fluctuations of g x , g y , and g z . Unless mentioned 
explicitly, we assume that g x , g v , and g z share the same distribution. 

• The dipole-dipole coupling strength D = 51.88 GHz, which is fixed. 

• The concentration n of magnetic impurities on the diamond lattice. 



Appendix B: Numerical solution of the phenomenological model 



As in the case of the Bloch equations, if the relaxation time T\ is finite, it is useful to be able to specify both the 
initial value (S(t = 0)) = (S(0)) of the magnetization and its stationary-state value (S(t = 00)) = (S)o- Therefore, 
we extend Eq. (gHJ) to 

|<S(*)) = A(S(i))+b, (Bl) 

where 

/ -1/T a 27rf*f \ 

A = -2ir£ z F -I/T2 irh p (l + ()(2 + +£y)F R , (B2) 

V -7r/ lp (i + o(2 + e + e y )^R -1/T1 J 

and b T = (S) /Ti. The formal solution of Eq. (|BH reads 

(S(t))(t) = e tA (S(0))+ / e uA hdu 

Jo 

= e tA (S(0))+A- 1 (l-e tA )b. (B3) 
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We integrate Eq. (|B1|) . that is we compute e tA , using the product-formula e tA — ( e rA i/ 2 e rA 2 e rA 1 /2\ m + C(r 3 ) (3C 
where r = i/m, A = Ai + A2 and 

-1/Ta 
Ai = I -1/T 2 

-l/Ii, 

2tt£ 2 F 

A 2 =1-2^0 ^(1 + 0(2 + ^ + ^)^1- (B4) 

-^ p (l + C)(2 + ?"+e)F R 



In detail, we have 



e x Al 



'e- 1 /^ 
e" 1 /^ 
e- 1 ^ 1 , 



1 - (6/fi) 2 (l -cosrO) (6/n)sinrfi (a6/^ 2 )(l - cosrO) 
rAa = ( -(6/n)sinrf2 cosrft (a/rj)sinrft j . (B5) 

(a&/^ 2 )(l -cosrfi) -(a/Q)sinrfi 1 - (a/fl) 2 (l -cos-rfi) 



where a = 2tt£ z .Fo, & = nh p (l + 0(1 + £ x + £ y )^R, and 17 = (a 2 + & 2 ) 1 / 2 . 
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